%% it maps the ISS subsidy to parameters of the model

y20      =19400
bISS     =60000
zeta_n   =70000
bnk_rate0=0.075;

k20      = y20*k_s_ben/(y_s/l_s)
b20      = k20-varzeta_s/(delta+r)/k_s_ben*k20

xb       = (coupon_s+1/8)/(r_c+delta_s+bnk_rate0+1/8)
xg       = (coupon_s+delta_s+bnk_rate0+1/8)/(r_c+delta_s+bnk_rate0+1/8)
zeta_g   = (xg-xb)*bISS/xg
zeta_i   = (coupon_s)/(r+delta_s+bnk_rate0+1/8)*bISS/xg
sub_size = (zeta_g+zeta_i+zeta_n)/k20
%sub_time = (zeta_n/2+zeta_i)/(zeta_g+zeta_i+zeta_n)
sub_time = (zeta_n+zeta_i)/(zeta_g+zeta_i+zeta_n)

